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ABSTRACT 

We study the soft X-ray variability of Cygnus X-3. By combining data from the All- 
Sky Monitor and Proportional Counter Array instruments on the RXTE satellite 
with EXOSAT Medium Energy (ME) detector observations, we are able to analyse 
the power density spectrum (PDS) of the source from — 0.1 Hz, thus covering 
time-scales from seconds to years. As the data on the longer time-scales are unevenly 
sampled, we combine traditional power spectral techniques with simulations to analyse 
the variability in this range. The PDS at higher frequencies (> 10 -3 Hz) are for the 
first time compared for all states of this source. We find that it is for all states well- 
described by a power-law, with index ~ — 2 in the soft states and a tendency for a less 
steep power-law in the hard state. At longer time-scales, we study the effect of the state 
transitions on the PDS, and find that the variability below ~ 10~ 7 Hz is dominated by 
the transitions. Furthermore, we find no correlation between the length of a high/soft 
state episode and the time since the previous high/soft state. On intermediate time- 
scales we find evidence for a break in the PDS at time-scales of the order of the orbital 
period. This may be interpreted as evidence for the existence of a tidal resonance 
in the accretion disc around the compact object, and constraining the mass ratio to 
M 2 /Mi < 0.3. 
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1 INTRODUCTION 

The persistent X-ray source Cygnus X-3 was discovered 40 
years ago and is one of brightest and most extensively stud- 
ied X-ray binary (XRB) systems. It is however a peculiar 
source that still eludes simple classificatio n. It is a close bi- 
nary with an orbital period of 4.8 hours (jParsignault et all 
1 19721). The identification o f the donor as a Wolf-Rayet star 
l|van Kerkwiik et alj Il992f) classifies it as a high-mass X- 
ray binary (HMXB). The strong wind from the Wolf-Rayet 
companion, surrounding the entire system makes radial ve- 
locity measurements difficult and the masses of the com- 
ponents are therefore uncertain. Recent studies based on 
the luminosity at the observed state tr ansitions suggest that 
the accreto r is a ~ 30 Mq black hole l|Hialmarsdotter et al.l 
120081 . 120091 ). This would make Cyg X-3 similar to the two 
ot her Wolf-Rayet X-ray binaries discovered so far, IC 1 X- 
1 (|Prestwich et al.ll2007l; ISilverman fc Filippenk3l2008h and 



NGC 300 X-l jCarpano et al.l l2007al lbh . where at least one 
of them, IC 10 X-l, contains a black hole of mass 23-34M0. 
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1.1 Spectral states 



Cyg X-3 displays five distinct spectral states 
dSzostek fc Zdziarskil l2004h . similar to those of the 
XRB transients (Cyg X-3 itself is however a persistent 
source), but with the hardest state peaking at much lower 
energies. The presence of strong absorption has made 
it difficult to interpret the intrinsic spectral shapes and 
whether they represent true state tra nsitions or are merely 
an eff ect of variable local absorption. IHialmarsdotter et al.l 

(2008) , however, presented several strong arguments that 
the transitions indeed correspond t o an intrinsic change 
in th e accretion geometry, and in IHialmarsdotter et al.l 

(2009) RXTE data from all spectral states were consis- 
te ntly modelled. In this pape r we use the nomenclature 
of IHialmarsdotter et al.l (|2009l ) to describe the states as 
the hard, intermediate, very h igh, soft non-therma l and 
ultrasoft (see figs. 1 and 7 in IHialmarsdotter et al.l 120091 
for spectral shapes and envisaged accretion geometries) in 
order of increasing softness. The very high, soft non-thermal 
and ultrasoft states will for most of the time be combined 
and referred to as the "soft state" or states as they differ 
mainly in the shape and strength of the high-energy tail. 
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1.2 Variability 

The X-ray temporal properties of Cyg X-3 are not as well 
studied as those of for example Cyg X-l. The strong orbital 
modulation with a peri od of 4.8h was discovered in 1972 
l|Parsignault et al]|l972h . In a study of the short time vari - 
ability using EXOSAT data, Ivan der Klis fc Jansenl (|l985T ) 
reported transient quasi-periodic features in the lightcurves 
during the soft states, and possibly also in the hard state. 
The periods for these features were in the 50-1500 s range, 
with a given feature s urviving fo r 5-40 periods. Similar re- 
sults were reported bv lRao et all l|l99ll ). who found a tran- 
sient oscillation with a ~ 121 s period. 

A po wer density spectrum (P DS) of Cyg X-3 was pre- 
sented by Willingal e et al.l Jl985), covering the frequency 
ranee ~ 10" 5 -0.1Hz. It was found that it was well described 
by a power-law of index —1.8, with a flattening at lower 
frequencies. A strong component was seen from the orbital 
modulation. No variability signal was detected above 0.1 Hz. 
A more detailed study of the X-ray variability above 1 Hz by 
iBerger fc van der Klisl (|l994h found an upper limit of 12% 
rms on the var iability in this ran ge. 

Recently, IChoudhurv et ail (|2004T ) presented results 
from analysing RXTE/PCA observations of Cyg X-3, in- 
cluding corrections for the binary variations. They compare 
observations made at three different times, and find that the 
PDS does not change significantly; above 10 -4 Hz it can be 
fitted by a power-law with index —1.5. Using RXTEj ASM 
data, they find that the power-spectrum on longer time- 
scales (> 1 day) is flatter, and through extrapolation predict 
a break around 10 -4 Hz. 

In this paper we study the PDS of Cyg X-3 in the 
range 10" 9 - 0.1 Hz using RXTE/PCA, RXTE/ ASM 
and EXOSAT data. Based on the spectral analysis of 
iHialmarsdotter et alJ (|2009h we look for systematic differ- 
ences in the PDS between the different states. We also study 
the long-term variability and the full broad-band PDS be- 
tween the shortest and longest time-scales available. The 
data and methods used are described in Sect. 2, and the re- 
sults are presented in Sect. 3. This is followed by a discussion 
in Sect. 4, and in Sect. 5 we summarise our conclusions. 



2 DATA AND ANALYSIS 



Levine et al.l 
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instruments onboard the RXTE satel- 



lite. The RXTEj ASM provides continual coverage of the 
X-ray sky in the 2-12 keV range. A given source will be 
observed several times each day, and the data is available 
either in dwell-by-dwell format, or one-day averages. The 
RXTE/PCA is used for pointed observations, typically last- 
ing up to about an hour. The pointing data i n this pa- 
per is the same as in IHialmarsdotter et alJ (2009), and that 
spectral analysis is used when determining the state of the 
source. 

To create the PDS, PCA lightcurves were extracted 
with 100 ms resolution in the 2-9 keV band. The lightcurves 
were then separated in intervals of 4096 bins, and a PDS was 
computed for each interval. Only intervals without gaps were 
retained, and the resulting PDS were then averaged to form 



the final power spectrum. Since the data stem from many 
different observations with varying setup, the Proportional 
Counter Units (PCUs) in use vary between observations. All 
lightcurves have therefore been normalized to one PCU. 

To bridge the gap in frequency coverage between the 
PCA and ASM datasets, we have also included data from the 
EXO SAT /Medium Energy (ME) instrument (|Turner et all 
1981). The data comprise the 19 observations of Cyg X-3 
made with EXOSAT spanning more than 9ks (in total 22 
obs ervations were made of Cyg X -3 between 1983 and 1985, 
see IBerger fc van der Klisl fl 994 for a log). Due to gaps in 
the data, these observations were divided into 26 segments, 
spanning between 3000 and 27 750 s. The data cover the en- 
ergy range 0.8-8.9 keV. Include d in this dataset i s the l ong 
(106 ks) observation analysed by IWillingale et al.l l|l985h . 

X-ray lightcurves of Cyg X-3 show strong modulation at 
the 4.8 hour orbital period. We have corrected all PCA, ASM 
dwell-by-dwell and EXOSAT lightcurves for this modulation 
by dividing fluxes with a modulation template, taken from 
van der Klis fc Bonnet-Bidaudl (ll989T ). The amplitude of the 
template wa s fitted to the data with phases give n by the 
ephemeris in Ivan der Klis fc Bonnet^ Bidaudl (|l989h . 

A problem with interpreting ASM data is dealing with 
systematic effects. The data show variations with periodic 
components on time-scales matching one year, one day and 
one orbit. These periods are related both to variations in 
the observing conditions and in the instrument properties. 
Although correction factors are applied when producing the 
ASM lightcurve data, it is not possible to completely elim- 
inate these effects. As a result, some systematic effects will 
still remain; however, for a bright source such as Cyg X-3, 
these systematics are only present on a low level for long 
time-scales (< 3%; A. Levine 2008, priv. comm.). 



2.1 Unevenly sampled data 

Although the observations made by the RXTE/ ASM are 
combined into one-day averages, they are not evenly sam- 
pled. A given source will be observed several times per day, 
in principle allowing variability on shorter time-scales to be 
detected by analysing the dwell-by-dwell data. In the case of 
Cyg X-3 it also allows us to correct the data for the orbital 
modulation. 

The RXTE/ 'ASM data was split up into different states, 
based on the flux level. Small pieces of data in-between the 
defined states were not used. The data for each state was 
then divided into segments with lengths between 250 and 
1000 points. Power density spectra were computed with a 
Lomb-Scargle algorithm, normalized to variance per inten- 
sity squared and frequency unit. A logarithmically averaged 
PDS was computed and a power-law fitted to the frequency 
range 2 x 10" 7 -6 x 10" 6 Hz. 

The systematic effects on the PDS caused by the data 
sampling were investigated by Monte Carlo simulations. 
Synthetic data sets, with the same time sampling as the 
observations were produced and analysed in the same way. 
The simulated lightcurves were produced by shot noise. By 
injecting exponential shots with different distributions in de- 
cay times we simulated lightcurves for which the PDS were 
power- laws with index ranging from —1 to —2. The differ- 
ence between the simulated and measured power-law index 
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Figure 1. Power spectra from RXTE/PCA data for all states of 
Cyg X-3, shifted for clarity. The dashed line is a power-law with 
index —2. The flattening at high frequencies is due to the Poisson 
level not having been subtracted. 
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Figure 2. Example fit of the PDS of Cyg X-3 in the soft non- 
thermal state (solid line). Rather than subtracting the Poisson 
level, we used a model comprising a constant component (dotted 
line) and a power-law (dashed line) to fit the data. 

allowed us to estimate the systematic errors introduced by 
the data sampling. 



3 RESULTS 

3.1 PDS at high frequencies: 10 3 0.1 Hz 

Using RXTE/PCA data, we have conducted a systematic 
study of the high frequency PDS of Cyg X-3. Figure [1] shows 
example PDS for all states seen in the source as defined by 
iHialmarsdotter et al.l (|2009l ). The PDS is in all cases well 
described by a steep power-law up to ~ 0.1 Hz. At higher 
frequencies, the signal is too weak to be detected. 

In order to test for systematic differences, we fit each 
PDS to determine the power-law index. Rather than sub- 
tract the Poisson level, we include it in the model as a con- 
stant. An example fit is shown in Fig. [2] 

The resulting indices are then grouped according to 
state. This allows us to compare the variations within each 
state to that between states. Table [1] lists results and uncer- 
tainties for each state, obtained by combining all individual 
fits for that state. As the number of observations in each 
state varies, the rightmost column gives the number of PDS 
combined in each state. We see that the index in all states is 



Table 1. Fitted power-law indices in the PDS of all states. The 
data are from RXTE/PCA, and the frequency range is 10" 3 -0.1 
Hz. 



close to —2. However, there appears to be a tendency for the 
power-law in the hard state PDS to be slightly less steep. In 
this state we also noted that there were larger fluctuations 
between the individual observations. The results are consis- 
tent with the index being the same in all the soft states. 

A power-law with index —2, sometimes referred to as 
red noise, can result from a random walk process. Our re- 
sults could thus indicate that such a process dominates the 
high-frequency power spectrum in all states of the source, 
drowning any intrinsic variability. This will be further ad- 
dressed in Sect. 14.21 

3.2 Long time-scales: 10 9 10 5 Hz 

The RXTE/ ASM has observed Cyg X-3 for over 12 years, 
giving unprecedented opportunity to study the long-term 
behaviour of the source. Figure [3] shows the 2-12 keV 
lightcurve of Cyg X-3. The state transitions between the 
hard state with low ASM flux and the soft states with high 
ASM flux are clearly visible. 

We now calculate the PDS using the method described 
in Sect. 12.11 To check for differences between the states, we 
first split the data into segments based on the RXTE/ ASM 
count level: hard state (count rate ~ 10 cts/s) and soft state 
(~ 20 cts/s). A PDS was calculated for each segment, and 
all PDS of the same state were then averaged for the final 
result. As for the higher frequencies, the PDS in both states 
is well described by a power-law. The index in the hard 
state is -0.93 ± 0.09 and for the soft state -0.78 ± 0.07. 
However, these PDS are distorted by the effects of the data 
sampling, and the indices derived do not give a true rep- 
resentation of the intrinsic variability. From our approach 
using simulations, we find that the power-law of the intrin- 
sic PDS is steeper: the observed results correspond to an 
intrinsic power-law index of —1.3 in the hard state and —1.5 
in the soft state. The relative error of the corrected indices 
are roughly the same as for the uncorrected ones, i.e., about 
10%. The difference in index between the two states is there- 
fore not significant. 

3.2.1 State transitions 

We have also combined all data, and calculated the PDS for 
the entire ASM lightcurve. The resulting index is in this case 
steeper than when either of the two states is treated indi- 
vidually. Looking at the RXTE/ ASM lightcurve, it is clear 
that the state transitions are responsible for much variabil- 
ity on long time-scales. To determine the contribution of the 
transitions to the PDS we produced an artificial lightcurve 
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Figure 3. 2—12 keV lightcurve of Cyg X-3 . The data span more than 12 years, from March 1996 to July 2008. The soft states are 
clearly visible as periods of increased flux. 
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Figure 4. Results of testing the contribution of state transitions to the PDS of Cyg X-3. Left panel: The artificial lightcurve created 
with only three levels, corresponding to low, intermediate and high states. Right panel: The PDS of the artificial data (black curve) 
overlaid on that of the real lightcurve (grey curve). The two curves agree almost completely below ~ 10~ 7 Hz. 



with only three levels, corresponding to the hard, interme- 
diate and soft states. This removes all variability within 
each state, isolating the power produced by the transitions 
themselves. The resulting lightcurve is shown in Fig. [4] (left 
panel). The artificial time series was evenly sampled, with 
a bin length of 0.1 days and with no data gaps. We then 
computed the PDS for our artificial data and compared to 
the real PDS. The two curves are shown overlaid in the right 
panel of Fig. 0] 

In the figure we see that there is almost complete agree- 
ment between the real and artificial data at lower frequen- 
cies. We stress that we have not performed any fitting of the 
two curves. Nevertheless, they follow each other almost com- 
pletely up to ~ 1CT 7 Hz. The PDS of Cyg X-3 below this 
frequency can therefore be explained simply by the variabil- 
ity of the ASM flux connected with the state transitions. 
It is also clear that the time sampling of the ASM data is 
not responsible for the structure of the PDS at low frequen- 
cies since these effects are not present in our artificial state 
lightcurve. 

Using the ASM data we have looked for possible corre- 
lations between the time that Cyg X-3 stays in a soft state 
and the length of the preceding and of the following hard 
and intermediate state. With the 15 soft states available for 
this analysis no obvious deviation from a random scatter 
was found. 

We also analysed the distribution of the length of time 
that the source remained in a hard or soft state. The mean 
lengths of the 19 hard and 15 soft states were 103 and 
98 days, respectively. The Cumulative distribution function 



for each of the states was compared with that of the wait- 
ing times for a Poisson process with the same mean value. 
The Kolmogorov-Smirnov D-value was 0.125 and 0.131 with 
corresponsing significances of 0.907 and 0.943, respectively. 
Both distributions are therefore consistent with a process 
where the probability of a transition at any time is indepen- 
dent of the time since the last transition. The fact that the 
distributions are well explained by a Poisson process does 
not exclude other possibilities. 

3.3 The complete picture 

To complement the data from the RXTE/ ASM and 
RXTE/PCA instruments, we use the data from EXOS AT tor 
the intermediate time-scales (10~ 5 — 10 -2 Hz). As expected, 
the PDS at intermediate frequencies is well described by a 
broken power-law: at lower frequencies the index matches 
that of the ASM data, and at higher frequencies the power- 
law steepens, matching the PCA data. 

Of the 19 observations used, 3 are in the hard state, 1 
is in the intermediate state and 15 are in the different soft 
states (which we again combine together into one PDS). The 
resulting indices are presented in Table [2] We find a slightly 
higher index in the hard state than in the soft and interme- 
diate states. Although the difference is not significant, this 
matches the trend seen in the PCA data, where the hard 
state shows a slightly less steep power-law. 

In all our EXOSAT observations we find signs of a break 
around 10 -4 Hz. The break is seen in the PDS of both the 
uncorrected lightcurves and of the lightcurves with orbital 



© 0000 RAS, MNRAS 000, 000-000 



X-ray variability of Cygnus X-3 5 




10-* 

Frequency [Hz] 



Figure 5. Broadband power spectrum of Cyg X-3. The solid 
line is a power-law with index —2.1, fit to the RXTE/PCA data, 
and the dashed line is a power-law with index —1.3 matching the 
RXTE/ ASM data. The dotted line indicates the Poisson level. 



State 


Index 


Hard 


-1.78 ±0.07 


Intermediate 


-1.9 ±0.1 


Soft 


-2.00 ±0.04 



Table 2. Fitted power-law indices above 10 4 Hz in three states 
using EXOSAT data,. 



demodulation applied. The limited frequency range below 
the break in the EXOSAT data together with the orbital 
modulation and demodulation make it difficult to determine 
the low-frequency index. However, t he long (110 ks) interme - 
diate state observation presented in I Willingale et all (| 19851 ) 
does allow a fit at lower frequencies. Analysing this observa- 
tion alone, we find an index of -1.08 ± 0.53 below 10" 4 Hz. 
From our simulations we can determine that the break it- 
self not significant; however, extrapolating the PDS at low 
and high frequencies gives an intersection around 10 -4 Hz. 
Fig. [S] shows the total combined broad-band PDS covering 
the frequency range 10" 7 - 0.1 Hz. The RXTE/ ASM and 
EXOSAT data combine all soft states, and the RXTE/PCA 
data shown is from the soft non-thermal state. 



4 DISCUSSION 

4.1 The PDS as a state indicator 

In many XRB sources, the power spectrum at high frequen- 
cies changes with the sp ectral state of the source. As shown 
by several authors (e.g ., iBelloni et al. 1 120021 ; lAxelsson et al.l 
l200rj ; lDone et al.|[2007n studies of the PDS complement spec- 
tral analysis, and may in some cases provide a clearer indi- 
cation of state than the spectral energy distribution. Unlike 
other X-ray binary sources, for exmple Cyg X-l, the PDS 
of Cyg X-3 do not show large variations between different 
spectral states. The PDS in the different soft states are indis- 
tinguishable from each other and the only possible difference 
we may see is a slight tendency for a less steep power-law in 
the hard state. 

The case for a true state transition in Cyg X-3 would 



be considerably strengthened had quasi-periodic oscillations 
(QPOs) been detected in the hard state PDS, as in other 
sources, but we see no evid ence of such features (nor of any 
transient oscillations, as in Ivan der Klis fc Jansenll 19851 ). In 
the high-mass black hole binary Cyg X-l, these features 
dominate the hard state PDS above ~ 0.01 Hz, with the 
lowest frequency feature peaking at 0.1-1 Hz. In this range 
the PDS of Cyg X-3 is however below the Poisson level which 
would prevent such features from being detected. Although 
the origin of the QPO features is not yet understood, many 
models predict an inverse scaling with mass (i .e., a scaling 
with the Keplerian frequency, see for example Ivan der Klisl 
l2004h . 

A generally accepted model for the accretion flow 
around the compact object is a two-phase flow, with an 
outer cool thin disc that is truncated at some radius and 
replaced by a hot inner flow. The truncation radius is deter- 
mined by the mass accretion rate - at high accretion rates 
the thin disc may extend in to the inne rmost stable orb it 
(for a thorough review of the model, see iDone et al.ll2007T ). 
This model provides a natural candidate for a specific ra- 
dius where the frequencies are realised: the truncation ra- 
dius. The low-frequency feature can be estimated using the 
viscous time-sc ale (the longest t ime-scale) at the truncation 
radius, giving l|Done et al.ll2007l ): 



(i) 



where r is the truncation radius in units of gravita- 
tional radii and m the bla ck hole mass in solar masses. 
iHialmarsdotter et al.l (|2009l ) concluded that for the transi- 
tion to a hard state in Cyg X-3 to occur at the same Ed- 
dington luminosity as in other non-transient sources, the 
mass of the compact object must be ~ 30 Mq, three times 
that of Cyg X-l. Scaling up the low frequency component 
in the PDS of Cyg X-l for this mass , we would expect this 
component to appear at ~ 0.03 — 0.3 Hz. In this range, it 
could be visible if it is there. However, the exact truncation 
radius and blackbody temperature in Cyg X-3 are poorly 
known due to both absorption and emission features in the 
wind that affects the spec trum strongly at low energies (e.g. 
ISzostek fc Zdziarskill2008 | ) . Ac cording to the best-fit models 
of IHialmarsdotter et al.l (2009) the disc in the hard state of 
Cyg X-3 is truncated further in than in the hard state of 
Cyg X-l. It is therefore likely that the component would 
show up in the higher end of the frequency range, which is 
again close to where the variability drops below the Poisson 
level. In conclusion, the non-detection of any QPO features 
in the hard state PDS of Cyg X-3 neither confirms nor ar- 
gues against a true state transition in the source. 



4.2 Effects of the wind on the PDS 

Continuing the comparison with other XRBs, it is clear 
that the power-law in Cyg X-3 is quite steep - the most 
common value in other sources is an index close to —1 (at 
least in the soft st ates where the power-law dominates; see 
Ivan der Klisl |2004| . and references therein) . An interesting 
question is therefore what process could give rise to this 
PDS. As the system is enshrouded in heavy stellar wind, it 
is tempting to explain the PDS as the result of scattering 
in this medium. Such a scenario would predict that frequen- 
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cies above the random walk time-scale are suppressed in 
the PDS, much like the highe r frequencies in Cyg X-3. The 
explanation was suggested bv lBerger fc van der Klisl l| 19941 ) 
to explain the weak or missing high-frequency variability in 
Cyg X-3. If the wind optical depth does not vary much be- 
tween the states of the source, this could also explain the 
more or less constant appearance of the PDS. 

This scenario requires the Wolf-Rayet wind to be op- 
tically thick to scattering. The typical scattering time-scale 
for a spherical homogeneous sca ttering cloud of radius R 
and optical depth r is 0.5Rt/c (|Kvlafis fc Phinnevlll989l) . 
The ~ —2 index seen in the PDS of the PCA data con- 
tinues to 10 -4 Hz. Associating the break at 10 -4 Hz with 
the scattering time-scale and assuming a radius of the dense 
wind region of 10 12 cm (an estimate based on IR observa- 
tions; iFender et al] [l999), would however require an unreal- 
istic or>tiral_dfip_H2_of^600 ; 

ISzostek fc Zdziarskil (120081 ) modelled the Wolf-Rayet 
wind including absorption and emission features. They find 
an average r ~ 0.3 — 0.4. This value gives a scattering 
time-scale of the order of only a few seconds for a radius 
of 10 12 cm, or requir es a radius of 3 x 10 4 R0 to explain the 
break in the PDS. In lHialmarsdotter et all (j2009) it was sug- 
gested that the low electron temperature in the hard state 
of Cyg X-3 may be a result of Compton down-scattering in 
an optically thick wind, either from the Wolf-Rayet star or 
in the form of a disc wind/outflow (this effect was however 
not included in the modelling of the X-ray spectrum). Such 
a wind would have to have an optical depth of r ~3-5 to 
down-scatter typical 100 keV photons to the observed 20 
keV peak. However even such a dense wind would still have 
to be extremely large, ~ 2500 Rq. 

Thus, while scattering in a nearby surrounding medium 
may explain the la ck of high-frequency variabi lity (> 0.1 
Hz) as suggested bv lBerger fc van der Klisl l| 19941 ). the ~ —2 
index of the PDS cannot be attributed to random walk in 
the Wolf-Rayet wind or any local scattering medium. 

Allowing for a very large size of the scattering cloud 
does not solve the problem. Although a large scattering 
halo has been detected around Cyg X-3 and the time delay 
from sc attered X-rays used to determine the distance to the 
source (jPredehl et al.l l2000l ) . this is the result of scattering 
in interstellar dust. The optical depth in this medium is thus 
<C 1, and the halo is not a sign of the source being extended. 
However, an ext ended feature has been found around the 
Cyg X-3 system jHeindl et alj|2003h . It does not surround 
the entire source, and has been interpreted as a localized 
remnant of earlier outflows. Recently, the feature has been 
shown to display flux and phase correlat ions with Cyg X-3 , 
making such an interpretation less likely (|McColloughll2008h . 
A possibility is that the feature is the result of a small scat- 
tering cloud within 2 kpc of Cyg X-3, and can therefore not 
explain the high-frequency power spectrum. 

Our results thus rule out the steep power-law as being 
due to scattering effects in a Wolf-Rayet wind, disc wind, 
local or line of sight matter, with any reasonable physical 
parameters. 

4.3 Long time-scales and transitions 

The PDS of many XRBs are well described by a power- 
law at low frequencies. In a survey of low-mass XRBs, 



iGilfanov fc Arefievl (2005) found that most of the studied 
systems displayed a power spectrum with index close to 
— 1.3 at low frequencies. The physical origin of this vari- 
ability is however no t yet determined. A commonly adopted 
framework is that of iLvubarskiil |l997), where disturbances 
in the outer region of the disc propagate inwards. Faster 
variations thus get superposed on slower variations, and the 
emergent power spectrum will have the shape of a power- 
la w. A physica l mode l base d on this scenari o was p roposed 
by iKing et al.l (|2004T ). and iTitarchuk et al.1 (|2007h showed 
that a two-phase model (cool disc and hot flow) could suc- 
cessfully fit the broadband PDS of Cyg X-l and Cyg X-2. 
The models predict a break in the frequency spectrum at 
time-scales corresponding to the longest variations, and a 
flat (white noise) PDS below this. 

IGilfanov fc Arefievl l|2005h associate the longest time- 
scales in the accretion disc with the viscous time-scale at its 
outer boundary. We do not detect any break or flattening 
at the lowest end of our frequency range, that could be as- 
sociated with an outer disc radius. On long time-scales, it 
is clear that the variability introduced by the state transi- 
tions dominates the power spectrum, making such a feature 
hard to detect in the PDS of the complete ASM lightcurve 
(cf. Fig [4j if the disc is larger than half the orbital separa- 
tion. Since Cyg X-3 is a wind-accreting system, it is likely 
that the disc significantly smaller. A feature in the PDS con- 
nected to the outer edge of the disc could still be present at 
higher frequencies. This possibility will be discussed in the 
next section. 

Even if dominant at our lowest frequencies, it is clear 
that the power from the transitions alone does not explain 
the PDS at higher frequencies, showing that we are seeing in- 
trinsic variability from the X-ray emitting region. The PDS 
above 10 -7 Hz can be explained within the framework pre- 
sented above, with disturbances in the accretion disc prop- 
agating inwards and ultimately modulating the X-ray flux. 
If the steep power law at high frequencies is not due to 
any scattering effect, i t means that the PDS is intrinsically 
steep. In the model of ILvubarskiil J 19971 ) a power-law index 
of ~ — 1 arises if perturbations in the disc are similar at all 
radii. If the perturbations are greater at large radii, the re- 
sulting power spectral index is steeper. This can be coupled 
to the disc being thickened at large radii, as proposed for 
the su ggested accretion geometries in iHialmarsdotter et al.l 
|2009l ). 

Our results further show that the probability of a state 
transition in Cyg X-3 is not dependent on the amount of 
time since the previous transition. This is contrary to what 
one might expect for transient systems, where an outburst 
consumes a large fraction of the accretion disc. The time- 
scale of an outburst is therefore dependent on the amount 
of matter in the disc, which should correlate with the time 
spent in quiescence, when the disc is built up. The state tran- 
sitions in Cyg X-3 thereby seem more directly coupled to the 
mass loss rate of the companion star. This means that the 
viscous time-scale for a reconfiguration of the accretion flow 
as a response to a change in the accretion rate is shorter than 
the typical time-scale for variations in the mass-loss rate of 
the companion (or capture rate of the compact object). 
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4.4 The break in the power spectrum 

Although the break around 10~ 4 Hz seen in the EXOSAT 
PDS is not formally significant, the fact that it is consis- 
tent with the intersection of the power-laws from the (high- 
frequency) PCA and (low-frequency) ASM power spectra 
speaks in favour of the break being real. The s ame feature 
was seen in the results of IWillingale et al.l (|l985i ). who anal- 
ysed a single long EXOSAT observation. The orbital mod- 
ulation has been compensated for in both our analysis and 
that of IWillingale et all (|l985h . Nevertheless, the break is 
quite close to the time-scale of the orbital period which is 
slightly suspicious. We note that the demodulation in all 
cases had a negligible influence on the PDS below 10 Hz, 
so the main uncertainty is not the existence of a break, but 
its precise location. 

If the break is interpreted as the feature connected with 
the viscous time-scale at the outer edge of the disc, this 
would mean a very small accretion disc. The corresponding 
radius is approximately 3 x 10 -2 Rq, which is unreasonably 
small, and this interpretation can therefore be ruled out. 

It is interesting to further compar e our results to the 
analys is of low-mass XRBs presented bv lGilfanov fe Arefievl 
(2005). They found that all compact systems in their study 
showed a feature in the PDS close to the orbital frequency. 
These systems also have l ow mass ratios, q = M2/M1 < 0.3. 
iGilfanov fe Arefievl (|2005h argued that for low mass ratios, 
it is possible to excite a 3:1 resonance in the disc. The reso- 
nance occurs if the angular frequency of the orbital motion 
of a particle in the disc is commensurate with the angular 
frequency of the orbital motion of the secondary and will 
lead to asymmetry of the disc and a possible feature in the 
PDS. For systems with a higher mass ratio, the disc is tidally 
truncated before reaching the resonance radius. 

Tying the break at 10~ 4 Hz observed in Cyg X-3 with 
the 3:1 resonance indicates that the mass ratio in this sys- 
tem is low, q < 0.3. Although the mass of the compan- 
ion Wolf-Raye t star is poorly known, it is expected to be 
> 10 Mq (e.g.. [Crowther|[2007l ). giving a black hole mass of 
~ 30 Mq . This value is in ag reement with that suggested by 
iHialmarsdotter et alj <|2009h . 



5 CONCLUSIONS 

We have conducted a systematic analysis of the PDS of 
Cyg X-3 over the frequency range 10" 9 -1 Hz. For the first 
time, we compare the PDS in all states of the source. 

On short time-scales, we confirm previous results that 
the PDS is well described by a power-law with index close 
to —2. By comparing the variations within each state to 
that between states, we find that only the hard state index 
shows a tendency to deviate from this value. We do not find 
any evidence of QPO features in the PDS, but this neither 
confirms nor contradicts the presence of a truncated disc in 
the hard state. Our conclusion is that the PDS is not a good 
indicator of state in Cyg X-3. 

On longer time-scales we find that the variability is 
dominated by the state transitions: the PDS below ~ 10~ 7 
Hz can be completely explained by this pattern. We also 
find that there is no correlation between the probability for 
a state transition and the length of time spent in the state, 



suggesting a direct feedback between mass-capture rate and 
spectral state. 

Combining datasets, we were able to create a power 
spectrum covering the entire 10 _7 -0.1Hz range, and found 
evidence for a break in the PDS around 10 -4 Hz. This may 
be interpreted as evidence for the existence of a tidal reso- 
nance in the accretion disc around the compact object. 
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